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1. Introduction 



The possibility of arbitrary statistics continously interpolating between the Bose and the 
Fermi case in space dimensions less than three was suggested some years ago in a beautiful 
and pioneering paper [1] by Leinaas and Myrheim. A physical model of particles with 
arbitrary, or fractional, statistics was put forward by Wilczek [2], who baptized such objects 
as "anyons" . Anyon dynamics has lately received much attention (for comprehensive recent 
reviews and guides to the literature see ref. [3]- [5]), as it is thought to be a good candidate 
for the explanation of such important collective phenomena in condensed matter physics as 
the fractional quantum Hall effect and high temperature superconductivity. To this end, 
it is then crucial to investigate the dynamics of the many-anyon systems. 

The symmetry group of the multi-anyon state vectors is not the group of permuta- 
tions, it is the braid group. Due to this fact, even in the non-interacting case the many- 
anyon problem is not separable, and the exact solution has been obtained (by Leinaas and 
Myrheim in their original paper [1]) only for the case of N = 2 non-interacting anyons 
in a harmonic well. For the three-anyon problem Wu has obtained a special set of exact 
solutions [6] which have been recently generalized to the arbitrary iV-anyon system in a 
magnetic field by Chou [7] and by Dunne et at. [8]. In fact, the problem of N anyons 
moving in an external harmonic potential can be easily translated in that of iV charged 
anyons moving in a constant magnetic field. All the enrgy levels obtained from the special 
exact solutions [6]- [8] exhibit a strict linear dependence on the statistical parameter a, as 
in the exactly solved 2-anyon problem [1]. 

Recently, very important experiments have been performed by Sporre, Verbaarschot and 
Zahed, who have solved numerically the Schroedinger eigenvalue problem for the N = 3 
and the N = 4 anyon systems [9], [10]. The energy spectrum so obtained presents some 
levels with a slight non-linear dependence on a, and some with a strictly linear dependence 
on a. The latter include Wu's exact eigenvalues. 

In coincidence with the numerical investigations, Illuminati, Ravndal and Ruud have 
introduced a semi-classical approximation to determine the energy spectrum of the 3- 
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anyon system [11]. The energy levels are again linear in a and their slopes agree with the 
numerical ones. We have quite recently generalized the semi-classical quantization scheme 
to the general A-anyon system in external harmonic potential [12]. The semi-classical 
solution is obtained in connection with an exact limit of separability of the general many- 
anyon Hamiltonian, and it is shown to reproduce and explain all the main features of the 
numerical results also in the 4-anyon problem. 

In this paper we wish to address the many-anyon problem in external magnetic field 
in the framework of the exactly solvable model introduced in [12]. We will solve the 
eigenvalue problem from topological considerations specific of the classical trajectories of 
charged particles in a magnetic field. The energy spectrum so obtained reproduces not 
only the slopes (as in the case of the harmonic potential) but also the actual values of the 
numerical energy levels. 

2. Classical Trajectories in a Magnetic Field 

Consider the motion of a classical non-relativistic particle with mass m and charge e 
in a constant magnetic field of strength B. Assume the field directed along the ^-axis and 
the particle constrained to the xy-pl&ne (e.g. when the velocity along the field is zero). 
The particle then moves in cyclotron orbits with angular velocity uu = eB/m. Choosing 
the symmetric gauge A = |B x r, the lagrangian equation of motion for the particle has 
the general solution 



where r is the radius of the cyclotron orbit, R is the position vector of the center of the 
orbit, and i, j are the unit vectors along the x- and the y-axis respectively. The energy E 
and magnitude L of the angular momentum of the particle are 



r(t) = r [(coscut)i — (sinatf)j] + R, 



(1) 



E 




(2) 
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L = l -mu{R 2 -rl). (3) 

We see from eq.(3) that the classical orbits fall into two different topological classes, 
whether they enclose the origin of coordinates or not. When the angular momentum is 
positive the orbits do not enclose the origin of coordinates. Negative values of L correspond 
instead to orbits enclosing the origin. This feature of the classical orbits plays an important 
role in the semi-classical analysis of the many-anyon problem. 

3. The two-anyon system 

In the symmetric gauge, the Lagrangian for N anyons in a magnetic field is 

1 N 

L= 2 ^( m ^ + ef "i • B x r;) - ah ^ fa, (4) 

i=l i<j 

where = (xi,yi), and the azimuthal angle fa is defined by 

fa = arctan — — . (5) 

The statistical parameter a can be restricted to to take positive values in the interval 
[0, 1]. Bose satisfies is recovered for a = and Fermi statistics for a = 1. We could choose 
to reverse the sign of the topological term in the Lagrangian, but then a should vary 
in the interval [—1,0] to obtain the correct statistical influence of the winding numbers 
on the energy eigenvalue problem. The situation here is inverted respect to the case of 
the harmonic oscillator: there the topological term enters the Lagrangian with a positive 
relative sign, if a is chosen to be positive. 

We now specialize to the case N = 2. Introducing polar coordinates (p, 0), the separa- 
tion of the center of mass motion is immediate, and the effect of statistics is all contained 
in the Lagrangian for the relative motion 

77? • P 

L rel = -{f? + p^) + {-Bp 2 -ah)<p. (6) 
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The corresponding Hamiltonian is 



H = + fa + ah)2 - eB{n + ah) + -/ B2p2 ^ ^ 

Apart from the terms containing a, eq. (7) is exactly the Hamiltonian for an ordinary 
particle in a magnetic field. The eigenvalue problem is readily solved, in the same way as 
for the two-dimensional oscillator with statistical interaction originally treated in [1]. In 
terms of the Laguerre polynomials the eigenfunctions of H are 

ij n , m = e -^ 2 4 m+a| (e 2 )e im ^ |m+Q| , (8) 
where £ = p^Jmui/2%. The corresponding energy eigenvalues (modified Landau levels) are 

E n , m = (2n + 1 + \m + a\ - m - «)^; (9) 

the radial quantum number n is a positive integer, while the angular quantum number 
must be an even integer (positive or negative) if we model the two anyons as bosons with 
statistical interaction, and thus require the wave functions to be symmetric. The spectrum 
is drawn in Fig.l. From eq. (9) it follows that the energy is independent of the angular 
momentum and of a when m > 0. From the discussion of section (2) it corresponds to 
the classical relative orbit not encircling the origin of coordinates, i.e. to the two absolute 
classical orbits being separated and not winding around each other (see Fig. 2a). In turn, 
a zero winding number implies the corresponding quantum result that the energy does not 
depend on a. A second class of levels depending linearly on a with slope —2 is obtained 
when m < 0. The corresponding classical relative orbit encircles the origin, i.e. the orbit 
of one particle winds around the orbit of the other (Fig.2b), generating a non-zero winding 
number which, in turn, makes the corresponding energy eigenvlues depend on a. 

4. Separation of the many-anyon problem 

A semi-classical formula for the energy spectrum of N harmonically bound anyons was 
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derived in [12]. There the semi-classical result emerges as the exact solution of the quantum 
problem in a suitable limit of separability of the many-body Hamiltonian; it is the limit 
of non-crossing winding orbits. The same argument applies to the iV-anyon system in a 
magnetic field. Here we are motivated by the fact that the exact solution of the two-anyon 
problem corresponds to classically separated cyclotron orbits. We then assume that this 
property holds in general for N > 2. Introducing N — 1 Jacobi relative coordinates { p k } 
in polar form 



Pk = , (n + r 2 H hr fe _i - (k - l)r k ), 

, . V\ + V2 H \- Vk ~ ky k+ i 

(bu = arctan ■ 



x 1 + x 2 -\ h x k - kx k+ i ' 

k = 1,2,...,7V-1, (10) 

the separation of the N cyclotron orbits is expressed by the "clustering" of the Jacobi 
rays: p\ < p 2 < ■ ■ ■ < Pn~i- Here p\ is the relative ray of the two particles with the 
smallest separation, p 2 is the relative ray of the nearest particle to the first two, and so on. 
In the clustering approximation, comparing the definitions (5) and (10) of the azimuthal 
and the relative Jacobi angles, we have that 



J = 2,3,. ..,7V, (11) 

and the Lagrangian (4) for iV anyons expressed in Jacobi coordinates separates. The 
corresponding Hamiltonian for the relative motion reads 

1 N-l -, 1 

# = ^E (pI + -M, + M ) 2 - eB (p^ + kah ) + -/b 2 pD- (12) 

Z7Jl k=l Pk 4 

The eigenvalue problem is readily solved, and the multi-anyon wave functions ^/n factorize 
in the products of two-anyon wave functions 
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N-1 

^N= l[4>n k , mk , (13) 
k=l 

where 

i> nk ,m k = e -^^l7 +fe "l(4 2 )e imfc ^ei mfc+fca| , (14) 

with £ fc = pk\J muj /2h. The corresponding energy spectrum is the sum of two-anyon con- 
tributions with different statistical parameters 

N-1 

E = (2n k + 1 + \m k + ka\ — m k — ha) — . (15) 
k=i 2 

The radial quantum numbers n k are independent of each other and may take any 
positive integer value. The angular quantum numbers m k can take any positive or negative 
integer value, but they are however constrained to obey the quantum analogue of the 
clustering condition on the associated classical orbits. Since \m k \ gives the magnitude of 
the corresponding Bohr-Sommerfeld orbit p k , we get the relations 

\m k \ > |m fc _i|. (16) 

In addition, from the discussion of the two-anyon problem, mi can take only even integer 
values. 

Eqn.(15) combined with rule (16) gives the spectra drawn in Fig. 3 and Fig. 4 for the 3- 
and the 4-anyon systems, as a function of a. Neglecting the slight non-linearity of some of 
the numerical levels, we see that eqns.(15),(16) exactly reproduce both the slopes and the 
correct intercepts of the numerical eigenvalues [10]. A quick inspection of formula (15) tells 
that for a generic N, above a certain level there will be N(N — l)/2 levels with negative 
slope. 

The crucial ingredient in the solution of the many-anyon problem in a magnetic field 
appears to be the different topology of the classical orbits, encircling or non-encircling each 
other, which affects the angular momentum and thus the dynamics of the system. This 
is not true for the oscillator problem, where the separability of the many-anyon problem 

7 



corresponds only to classically non-crossing trajectories, with no distinction between en- 
circling and non-encircling orbits. The clustering rule (16) must then be relaxed in that 
case. Thus, for the oscillator problem the exact degeneracy is yet to be assessed. Work is 
in progress in that direction [13]. 

If we compare the eigenfunctions (14) with the special solutions obtained in refs. [6]-[8], 
we see that the latter are included as a subset of our solutions for extreme negative values 
of the angular momentum. Furthermore, it is easy to show in our semi-classical picture 
that the zeroes of the wave functions of Wu, Chou, Dunne, Lerda, Sciuto and Trugenberger 
correspond to non-crossing classical orbits. This explains qualitatively why we get them 
as a subset of our solutions. 

5. Discussion and Conclusions 

The approach to the many-anyon problem that we have put forward seems to have 
some strong justifications from the excellent agreement with the numerical results, from 
reproducing the exact solutions previously found, and also from the intuitive semi-classical 
picture. 

However, much work is to be done to assess the real connection between our solution 
and the complete solution of the general many-anyon problem. In particular, it is to be 
understood the role of the new wave functions that we have derived, which correspond to 
some of the linear and to all of the slightly curved numerical eigenvalues. This suggests the 
possibility of building the exact wave functions still missing through some perturbative or 
algebraic procedure applied to our solutions [14]. 

In conclusion, we wish to thank Pier Alberto Marchetti, Mario Tonin and Carlo Andrea 
Trugenberger for very stimulating discussions. 
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Figure Captions 

Fig.l. Energy spectrum of two anyons in a magnetic field, as a function of the statistical 
parameter. 

Fig. 2. Classical orbits of two particles in a magnetic field when the angular momentum is 
positive (2a), and when it is negative (2b). 

Fig. 3. Semi-classical energy spectrum of three anyons in a magnetic field. 
Fig. 4. Semi-classical energy spectrum of four anyons in a magnetic field. 
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